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TECHNIC VL MEMORANDUM 78164 


SIMPLIFIED MODEL OF STATISTICALLY STATIONARY 
SPACECRAFT ROTATION AND ASSOCIATED 
I NDUCED GRAV I TY ENV I RONMENTS 

I. INTRODUCTION 

Space vehicle rotations resulting from crew activity, thruster firings, 
etc. , and the associated induced g environments appear to be stochastic in 
character. Accordingly, to define an experiment to be performed aboard a 
space vehicle which is sensitive to vehicle rotations and the associated induced 
g environments, it would appear that statistical information concerning the risk 
of exceedance of critical rotation rates and g levels would be extremely useful. 
Clearly, the statistics of space vehicle motions and associated induced g environ- 
ments depend on the vehicle mass and geometry properties, the dynamic 
behavior of the crew (push-off, sneezing, etc.), the control system parameters, 
the mission profile, etc. Estimates of these statistics can be obtained from 
Monte Carlo simulation of the vehicle/forcing function system or from postflight 
analysis of g environment time histories acquired from onboard instrumentation 
(accelerometers, rate gyros) or look angle data. However, estimates of 
statistics of space vehicle motions and the associated induced g levels do not 
appear to be available at this time for the currently planned Space Transportation 
System (STS) missions with payloads involving g sensitive experiments (for 
example, those on Spacelab missions 1 and 3) other than estimates of typical 
and worst case rotations and associated g levels resulting from various kinds of 
discrete vehicle excitations. 1 Furthermore, statistical summaries of vehicle 
rotation and associated g environments measured on previous space flight 
missions do not appear available. However, a number of excellent reports are 
available on the effects of crew motion on spacecraft attitude. Reference 1 
documents the results of detailed simulations of the effects of crew motion on 
Apollo vehicle attitude, and Reference 2 provides the results of crew motion 
experiments conducted during the Skylab Program. This report attempts to 
provide estimates of exceedance statistics of vehicle rotations and associated 
g environments resulting from space vehicle motions. 

1. Lewis, R., Private Communication, 1978. 


The approach taken involves the use of basic assumptions concerning the 
statistics of the torque imposed on a spacecraft resulting from crew activity and 
thruster firings, i.c», that the imposed components of torque constitute a 
Gaussian process wherein the associated spectrum of torque in the frequency 
domain is that associated with a band-limited noise process with constant, 
nonzero spectral density over the frequency bandwidth of the process and zero 
spectral density for frequencies outside the bandwidth. The rigid-body equa- 
tions of motion are used to derive the statistics of the spacecraft rotations. 
Finally, the statistics of the associated g environments are derived from the 
vehicle rotation statistics by applying transformation formulae between inertial 
and rotating frames. The mathematical machinery of Rice’s theory of exceed- 
ances is used to obtain estimates of expected temporal rates of exceeding 
specified critical spacecraft rotation rates and associated g levels. Further- 
more, by assuming that the number of exceedances of rotation rate and induced 
gravity during an oi'bital experiment with duration time T arc Poisson processes, 
estimates of the risk associated with exceeding specified critical rotation rate 
and induced gravity levels at least once diming an experiment are obtained. It 
should be remembered that the calculations presented are speculative in nature 
and must await statistical analyses of results from Monte Carlo simulations 
and ''or of space vehicle acceleration data acquired from previous orbital 
missions. However, it is believed that the calculations presented are interesting 
and thought provoking and may bo useful to scientists and technologists who are 
developing space flight experiments which are sensitive to space vehicle 
accelerations. 


II. SPACECRAFT MOTION DESCRIPTION 

I.ct us consider a body-fixed frame of reference located at the spacecraft 
center of mass. This orthogonal frame of reference is fixed to the spacecraft 
with the x, y, and z axes directed along the principal axes of the spacecraft, 
bet us now consider a fluid container with center of mass located at position 
vector 7^ with components x^, y^, and relative to the spacecraft center of 

mass. Furthermore, assume the container is rigidly attached to the vehicle. 

It the vehicle undergoes rotation, then a fluid particle located at position 

with respect to the spacecraft center of mass with components x q , y Q , and z^ 

will experience a force per unit mass in response to the vehicle motion in ques- 
tion which is given by 


o 


» 


( 1 ) 



f = n x 7 { + jj x(nx 7 { ) - 2 nx u^r^t) +7 


where a is the inertial linear acceleration of the vehicle instantaneous center 
of mass resulting from a net force acting on the vehicle, iT is the rotation vector 
of the vehicle with components 0 , Q , and ft directed along the principal 

x y 2 

axes of the spacecraft, and (ft) denotes differentiation with respect to time [3]. 
The vector quantity u*(iy,t) is the velocity vector of the fluid particle with 

position vector r^ relative to a frame of reference with coordinates x', y\ and 

z' located at the fluid container center of mass and is fixed relative to the fluid 
container walls so that 


( 2 ) 


Figure 1 depicts the various frames of reference and the position vectors r^, 

r f , and r^’ . If the spatial extent of the fluid mass is small compared to the 

distance between the centers of mass of the spacecraft and the fluid container 
(i.e., Ir^’l << |7|), equation (1) may be written as 



x 

Figure 1. Principal axis and fluid container reference frames. 
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F = SI *7 + S2 * 7 ) + a* - 2s 2 V 7 ( 7 ’ t) 

o 'o' f 


( 3 ) 


The first three terms on the right side of equation ( 3) can be treated as a time- 
dependent gravitational body force per unit mass. The third term is a Coriolis 
force. Our analysis will be concerned with statistical definition of the gravita- 
tional-like body force terms in equation (3) (excluding the linear acceleration 
term) 


}f(t) s r + iT ' (IT V r ) 

' ’ o o 


(0 


and the angular velocity vector S2. It should be noted that if the fluid container 
is rotating (in addition to any rotations that may bo imparted by the spacecraft), 
then processional fluid flow analysis may be required, thus invalidating equa- 
tions (3) and ( I). This case is not treated in this report. 

The analysis of the statistics of git) with the linear acceleration term 
included increases the complexity of the analysis. The more general ease in 
which induced gravity results from both rotational and linear accelerations will 
lie analyzed in a suliscqucnt report. 

— ♦ — * 

The approximation r^, - r^ is extremely important because this approxi- 
mation moans that the gravitational body force which acts on a fluid particle duo 
to spacecraft rotations can be specified a priori in the sense that g to a 
sufficient degree of approximation is independent of the fluid flow dependent 
\ i aides. It r* cannot be approximated with r^, then "g cannot be specified 

■a : >*/,')•' because r ( . is a fiuid particle I.agrangian position vector which must 
then be determined as part of the fluid mechanics problem in question. 


A. Spacecraft Equations 

To develop the connection between the dynamics of the spacecraft and the 
vector quantities g and si , the rigid-body equations of motions tor the space- 
craft referenced to the principal axes are used. These equations, in vector 
notation, are given by 


-1 


( 5 ) 


i • ft = T+(n*i)xn 

where T is the net torque that acts on the spacecraft and 1 is the moment of 
inertia tensor with diagonal values 1,1 , and I and zeroes for the off-diagonal 
values 11] . x y z 


6. Linearized Equations 

The torque T results from rocket thruster firings and crew- activity. 
Some of the forces associated with crew activity (e.g. , push-offs, bending, 
instrument operation, etc.) have been quantified during Skylab missions and 
are documented in Reference 2, It appears that crew activity results in vehicle 
rotation with time scales on the order of seconds to a few tens of seconds. In our 
analysis, we shall be concerned with rocket thruster firing inputs required to 
keep a vehicle in a certain attitude in response to vehicle motions resulting 
from erew motions. These rocket firings produce relatively short period vari- 
ability in S? and hence in IT. Rocket firings associated with major changes in 
vehicle attitude will not be included in our analysis. Therefore, the terms 
involving products of rotation rates in equations (4) an I (5) may be neglected 
relative to the terms involving derivatives of rotation rates. It will be assumed 
that this is permissible and the following equations for ft and g will be used 


in the subsequent analysis: 


T . ft = T , 

(«) 

' = » " 7 o • 

(7) 


The validity of the decision to neglect the second-order terms involving ft to 
obtain these equations should be examined for each situation. The linearization 
process used to obtain equation (G) does not include the possible dependence of 
upon a control law and, in turn, a dependence of the controHaw upon ft . 

We shall bypass this issue by assuming that the statistics of T are known so 
that a statistical model of I? and can be developed via a rotational approach 
with equations (G) and (7). 





The use of equation (5) as a model to represent the combined effects of 
crew motion and rocket thrusters on spacecraft attitude is presumptuous. A 
model which includes crew motion exactly would be extremely complex because 
equations of motion for the vehicle and the crew members would be required. 

Those equations would include the effects of crew members attaching and 
detaching from the vehicle; vehicle crew member push-offs, sneezing, etc. , 
stochastic location of the crew members in time; and a host of other effects. 

An equation such as equation (5) would result from an analysis where IT would 
be the rotation rate of the vehicle about the vehicle center of mass (without crew) , 

I would be the moment of inertia tensor of the vehicle (without crew) , and T would 
contain the crew member/spacecrnft coupling terms and the torques imparted by 
the rocket thrusters. Equation (5) together with additional cciuations governing 
the crew would then require simultaneous solution for the dynamic dependent 
variables of the crew and spacecraft. Thus, equation (f>) should be viewed as 
an extremely simplified model of a complex system. However, it should be 
noted that equation (f>) is exact for the problem of calculating the response 
statistics of spacecraft motions resulting from stochastically imposed 
torques for a constant moment of inertia vehicle. 

It can be shown for spacecraft such as the Space Shuttle Orbiter with 
Spacclab as a payloadthat to neglect the terms in equations (1) and (G) which 
are second order in 52, we must require the lowest characteristic frequency 
ci’j. of the random process IT to be very large compared to any component of 

IT ; i.c. , 








C. Stochastic Models 

In this section we develop a stochastic model for Q and g based on 
equations (0) and (7) and as assumed stochastic model for T. 

1. Torque Stochastic Model . We hypothesize that the components of the 
torque vector arc mutually uncorrelated, statistically stationary Gaussian 
processes which have zero mean values and spectral density functions given by 


T ,T 

*T v ,T. (u) 2(7"-^”) • U 'l " ,UM ' U o 

, (9) 


,p (“') 0,0 

X* X 


< I u' I u' j 9 or < } u' I < 


with similar equations assumed for T and T . The quantity cr is radian 

y ^ 

frequency ( Paragraph Il.C. 2) , a is the standard deviation of T * and 

V x 

u’ and u> are upper :uui lower bound frequencies, respectively. The spocU’um 
o 1 

in equation (9) and all those that follow in the subsequent development are 
defined such that integration over the domain — w < u> < ■ v ' yields the auto- or 
cross-variance. The assumption that the components of the torque vector are 
uncorrelated appears to be reasonable for the crew activity contribution to 
torque. However, this assumption may not bo true for the thruster firing 
contribution to the torque vector. Nevertheless, this assumption is used in the 
analysis which follows, 

2. Rotation Vector Stochastic Process . We express il and T in terms 
of Fourier-Steltjes integrals HI 


H(t) = I e iut dH(u') 






~ ■ j 

i-iVj 



whore t is time and d„ (u ) denotes stochastic Fourier amplitude at frequency 
tv of the random process J\t) . Substituting equations (10) and (11) into 
equation ((>) yields the Fourier amplitudes of the V. process. For the \ 
component, we have 


i 



dS2 N (cc) 


dlMtc) 

iul 

\ 


(i-) 



i with similar equations for the y and / component rignW>ody equations Multi- 

i plying equation ( 12) by its complex conjugate evaluated at frequency u * and 

applying the ensemble average operator over all realizations of the (i2 f T) 
j process yields 




where the asterisk denotes complex conjugation and the angular brackets denote 
the ensemble average operator. The requirement of statistical orthogonaliU ot 
Fourier components or statistical stationarih of a random process in the time 
domain demands that 



where 

to 



is a spectral density function [ 1 1 • 


Thus, equation (11) reduces 



t 





*t T (w) 

X X 

“*I> 


(15) 


Combining equations (9) and (15) yields 


Vx (0,) 


= ST? 


T ,T 

x x 




, co < 1 0) I s to 

1 o 


o ( w ) ~ 0 » °< Iwl s> co, or to < |co| 

M tit 1 O 


X X 


(16) 


This result states that <f>^ ^ (co) cr oo“ 2 over the frequency domain in which 

_ (co) takes on nonzero values. Equations similar to equation ( 16) can be 
x x 

written for ft and ft . 

y z 


If it is assumed that the components of the torque vector are uncorrelated, 
it follows that the cross-spectral density functions of the components of U 
vanish so that 


i 


0 n ,« ^ ~ 0 ft ^ 
x y x y 




(co) = 0 
Z 


< CD < . (17) 


Statistical stationarity demands that the ensemble mean (ff } be equal to 
a constant vector IT Thus (ft} = 0 , which means <T>= 0 as hypothesized 

ah initio. For definiteness we shall set ft =0. 

o 

Integrating equation ( 16) over the domain - ® u ® yields the 
variance of ft^_ ; i.e. , 
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4 



« N ,« X 




i r J7 r 7 

\ o 


t 


( 1 «) 


where /> - iCj/w . Similar equations can be obtained for the remaining 
components of sT . 

In view of the linear relationship between 5,2 and T as expressed by 
equation ((») and the assumed Gaussian nature of the torque process, it follows 
that il is a Gaussian vector process. The previously mentioned data concerning 
the moments of the i 2 process provide sufficient information to determine any 
desired probability density functioned the components of S2 and hence to cal- 
culate any desired statistic of tho J2 process, 

3, Body Force Vector Stochastic Process . The Fourior-Stoltjos sto- 
chastic decomposition of g(t) is expressed mathematically as 





(ID) 


Substituting equations (|0) and (ID) into equation (7) yields 


tlg(u') di> ' r 

o 


C-'°) 


Application of tlio ensemble average operator to equation (20) (noting that 

<»> 0 and hence <"t> 0 ) yields 0 and hence ^g(t)^ 0. 

Thu;., tin* ensemble mean body force is zero in this model. 

Following the procedure outlined in Paragraph II. 0. 2, the following 
spectral density functions are obtained: 


*V s x 




( 21 ) 
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Vs M ' \<*j a) ’ 'V.‘ ! v, M 
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* \ z y y 
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Substituting the rotation rate* spectral density functions into equations (21) 
through (2(i) yields 


s x .s 
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with similar equations for g^ ami g^ and 
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with similar equations for the remaining cross-spectral densities and 
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Ileeause of the linear ^relationship between g and i.' 1 as expressed h\ 

equation (7)1 and because i! is a C mssian process, it follows that ?; is a 
(tans si an vector process. The previously mentioned data concerning the 





r .omenta of the g process provide sufficient infonnation to determing any 
• icsired probability density function of the components of g and hence to calcu- 
late any desired statistic of the g process. 

Since the g process is Gaussian, there exists an orthogonal transforma- 
tion of the spacecraft principal axis reference frame (x,y,z) to a new frame of 
r< ference (X, Y,Z) such that the cross correlations of the components of g 
i' the (X, Y,Z) frame vanish |5J. The components of g referenced to the 
(>:, Y,Z) frame are stochastically independent Gaussian processes. The trans- 
formation of the g process referenced in the (x,y,z) frame [i.e. , g(t;x,y,z)J 
to the (X, Y,Z) frame 11. e. , g(t;X,Y,Z)J is given by 

ff(t*X* Y, Z) = A * ^(t;x,y,z) , (35) 

where A is a 3 by 3 tensor with components that are functions of the auto- and 
cross-covariances of the g process referenced to the (x,y,z) frame, i.e., 
equations (28), (29), (30), (32), (33), and (34). The components of A can 
be determined by a relatively straightforward application of a three-way Euler 
angle transfort action [ 3] subject to the constraint that 


a 2 = a 2 = a 2 =0. (36) 

Sx»Ky e Z ,g X g Z' g Y 


Upon determination of the components of A , the spectral density functions of 
g (t;X, Y,Z) can cc obtained by Fourier transformation of equation (35) and 
formation oi 'ie appropriate square modulii. Thus, for example, Fourier trans- 
formation wf g^ yields 

dg = A dg (w) + A dg (w) + A dg (u>) . (37) 

x xx x xy y ’ xz z ' v ' 


Multiplying equation (37) by its complex conjugate at frequency w* and applying 
the condition of statistical orthogonality [equation (14) J yields 





A A X A \ V 


+ 2A A * (to) + A : o (. 0 ) 

xx xz g N ,H z v xy g y ,g y v ’ 

+ 2A A 0 (a*) + A* V ( 10 ) 

xy xz g ,g ' ' xz g ,g ' 1 

» / if 


Thus, the spectral density function of £ (t;X, Y.Z) c:m Ih' calculated directly 
from tl»c spectral density functions of ~g* (t;x,y,z) . The cross-spectral densin 
functions of g(t;X,Y, '/) vanish identically. The main point that we want to 
make is that there exists a coordinate system such that the components of g 
in that frame of reference are stochastically independent and that the statistics 
of g(t;X,Y,Z) are derivable from the statistics of g* (t;x,y, /.) . 

The spectral density functions of g“(t;X,Y,'/) are of the form 
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Similar expressions can be obtained for the Y and '/ components of g. 
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4. Cross Correlations of the Rotation and Body Force Vector Processes . 
The fact that g is derived from U implies a correlation between the g pro- 
cess with the If process. The following cross-spectral density functions can be 
calculated from equation (20): 
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Equations (41) through (43) show that parallel^components of dQ and dg are 
uncorrelated, so that parallel components of g and sj are uncorrelated. The 
cross-spectral densign functions of the mutually orthogonal components of and 
fl are complex according to equations ( 44 )_th rough (49). This means the 
mutually orthogonal components of g and fi are out of phase by ±90 deg. 

Thus, for example, the Fourier components of n lead those of g by 90 deg. 

y x 

Integrating equations (38) through (43) over the domain -® < to < 00 yields 


a 2 = <r i = 0 

V n y V 8 * 


(50) 


and likewise forjhe remaining cross correlations bejween the mutually orthogonal 
components of g and Q . These results show that g(t) is uncorrelated with 
Q (t) 1 The g and fi processes are, however, correlated for finite nonzero 
time delay t . This can be shown by the Fourier transformation of equations (44) 
through (49) to the t domain to obtain the cross-correlation function R(r) 
between the mutually orthogonal components of g(t + r) and i 2 (t). Thus, for 
example, 


R (t) ~ ~~ f e* aT <f> (a') d 

V n y 27r i \»« v 


u? 


Z rp CO T 

. o I t T o 

J— , iL_L r 

27rl } ( w -coj J * s ’ 

X ' O 1' CO jT 


(51) 


where the sine integral is tabulated in standard references. Thus, for r # 0, 
R (t) * 0 , with similar results for the remaining cross-correlation fune- 

fi tii 

X y _ _ 

tions between the orthogonal components of g :md S2 . 
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III. BODY FORCE PROCESS EXCEEDANCE STATI STI CS 

In this section we determine the exceedance statistics of each component 
of the g process referenced to the (X, Y,Z) frame discussed in Paragraph 
II.C.3. Our analysis will be concerned with a single component of g ; conse- 
quently, we will dispense with subscripts on g to denote components. We will 
use the symbol g to denote a component of If and a to denote the standard 
deviation of g . g 

A. Rice's Theorem and Expected Exceedance 
Rate of Body Force 

According to Rice (G,7J , for a stationary Gaussian process with zero 
mean, the expected number of exceedances of induced gravity per unit time 
which exceed level g is given by 


N 


g 


- gVZo* 

N e g 
g»o 


(52) 


The quantity N is the expected number of zero crossings of the quantity g 

from below and is related to the spectral density function 4> r (u.') through the 
following expression: g 


N = ~ 

g.o 2 7T 


-CO 


f a) 2 0g(a') da' 


V 2 


X 


(03) 


Thus, the substitution of equation (39) into equation (52) yields for our assumed 
stochastic g process 
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where 0= co^/to^ . The quantity fi serves as a relative measure ol the spectral 


bandwidth of the g jitter process. For a relatively broad-banded process in 

.s J 4, 


which B- 0, 2?rN /to = 3 /z ; and for a narrow process in which B 
g.o o 


1, we 


have 27rN /to =1. If the process were characterized by a monochromatic 
g.o o 

spectral density function, namely 


p(u>) 


g 


6(to 


- to ) 
o' 


+ 6 (to + to ) 
o' 


2 


( 55 ) 


where 6( ) is the Dirac delta function, then the zero-crossing rate would be 
given by = co q / 2 7r as calculated for [i - 1 with the g jitter spectral model 

given by equation (39). Thus, the model given by equation (39) includes broad 
banded and narrow-banded g jitter processes with the limiting case of mono- 
chromatic g jitter. 

Substituting equation (54) into equation (52) yields 


2,-N 

to 




v. 


-g 2 /2(T 2 

g 


(50) 


The nondimcnsional exceedance rate 2 ttN /to is plotted in Figure 2 as a func- 

g o 

tion of g/o for 0 and fi- 0.999. Thus, we conclude from Figure 2 that 
& 

for any given value of g/o^ the exceedance rate will vary by only a factor of 

1.733 over Uie admissible range of />. This means that for the model selected 

the exceedance rate N of the g process is only mildly dependent on the 

g 

parameter />, The dependence of N on />’ occurs through the zcx’o-crossing 

g 

rate N l equation ( 54) ] and, as px’cviously noted, 3 2 r- 2ttN /to s 1 

g.o • g,o o 
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Figure 2. Nondimcnsional expected rate 2 tN of exceeding the 

nondimcnsional level | g|_ ''it with positive slope for g > 0 or 

e g 

negative slope for g < 0, (The expected rate of exceeding 

level |g| o with positive or negative slope is 2N .) 
eg g 


for the full range of variation of the parameter p\ i.e. , 0 r- p 1. Thus, for 
example, if the g process has an upper-bound frequency of u' q = 2jt rad sec” 1 , 

then we have for the expected zero-crossing rate 0.577 < q 1 sec -1 . 


It should be noted that the g process can exceed critical values of |g| 

when g < 0 with negative slope; consequently, the rate at which the g process 

exceed a critical value is 2N . 

it 


f 



* 




] 
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B. Risk of Body Force Exceeding a Critical Value 


We now seek to determine the risk that the g process will exceed a i 

critical value for a given orbital experiment duration time. Clearly, the larger 

the duration time T of an experiment, the higher the risk a critical value of 

Igl (Igl ) wUl be exceeded. Ideally, we wish to know the probabilistic struc- 
c 

turc of the random time T when |g| is exceeded. This problem is called the 
’’first passage” problem. 

As previously noted, we consider the case where g is a stationary process. 

Furthermore, the g process is symmetrically distributed in the positive and 
negative ranges, and the upper and lower bounds arc also symmetrical. Thus, 
the probability of exceeding a critical g level, lgl c » at any given instant is 


Pl(g ^ Igl J n < K ? - Igl } 1 = 2P|g rr |g| J 
c c c 


(57) 


To estimate the risk associated with |g| > |gl c for a given experiment dura- 
tion time T , we shall make the arbitrary assumption that the exceedances of 

the |gi pi'oeess above level |g| arrive independently. We now denote by Q(T) 

c 

the number of exceedances of |g| at level |gl c over the experiment duration 

time T. Clearly, the process Q(T) is a Poisson process, and the probability 
of Q(T) being less than or equal to an assigned value (for example, q) accord- 
ing to Tin 18| , is given by 


P(Q(T) - q.T) c~ XT , (58) 


where A is a parameter. The probability of no exceedance of the g process 

above the critical value |g| in time interval T follows by setting q - 0 in 

c 

equation ( 58) , so that 


-XT 



P(Q(T) = 0,T) - e 


(59) 





By definition of the Poisson process wc set 


X =* 2N 

g 


Thus, 


-2N T 

P(Q(T) = 0,T) = e g 


((50) 


((51) 


Now, the risk R that the g process will exceed the critical value Igl^ at least 
once during an experiment of duration time T is 


-2N T 

R •- 1 - e p 


( 02 ) 


I 

f 

3 


Eliminating the expected exceedance rate 
yields 



between equations (50) and (02) 



(T 

S 



HKzA 
i - 7 


gfn (1 -R) 

Tu' 

o 


(03) 


This formula permits the calculation of critical g level (i.e. , lgl c ) » 8 a func- 
tion of risk R of the quantity Igl exceeding level (gl^ at least once during 

an experiment duration time T and the g environment spectral model parame- 
ters Li and u> . Figures 3 and 1 contain plots of Igl 'tr as a function of 
o eg 

Tu 1 for various values of risk for p 0 and 0.999, respectively, 
o 


It may be concluded from these figures that the variability in |g| due 

to variation in u' T and/or ii decreases as R decreases, 
o 
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C. Body Force Envelope Exceedance Rate 

The most questionable aspect of the previous analysis is the arbitrary 

assumption that the arrival of the threshold crossings of g above and below the 

critical levels |g| and -|g| with positive and negative slopes, respectively, 
c c 

are independent events. This assumption is especially unacceptable for narrow- 
band g jitter bee: use the threshold crossings of narrow-banded g jitter will tend 
to occur in clumps. Once there is a crossing of |g| over a threshold Igl » the 

probability is high that the tallowing excursion will produce another crossing. 
However, we note that the crossing of the same threshold by the envelope of 
the g process must precede the first crossing in each clump. Accordingly, when 
there are many excursions in each clump, the time of a threshold crossing by 
the envelope is nearly the same as the time of the first crossing in each clump. 
Thus, although it is more acceptable to treat the threshold crossings of an 
envelope of a narrow-banded random process as independent events, we can 
improve the analysis in the previous section by using the expected rate of 
threshold crossings of the envelope process for \ in equation (59). 

From Rice [6,7j, we assume the narrow-banded g jitter process can be 
expressed as 

g(t) = A(t) cos (u^t + 0(t)) , (64) 

where u is a representative wide-band frequency of the g process and A(t) 
m 

and o( t) are random processes which vary much more slowly than g(t) with 

respect to t. The process A(t) is nonnegative. Since the spectral density 

function of our assumed g jitter process is symmetric about the frequency 

a,' = cu ( 1 + li)/2 on the half interval 0 < w < ® , it is clear that « = w . 

so ms 

According to Rice, the random process A(t) is the envelope process of the 

g process. Since the g process is Gaussian, it can be shown that the expected 

exceedance rate or threshold crossing rate with positive slope of the envelope 

of the g process at level A is given by 


a a i i “A 2 /2cr 2 

M = — — T7 0 g , (65) 

°B <*r)* 






where 



j (u.'^ "a? 0 (u,') dit* 

mg' 

-iTl n 


(««) 


Substituting the g process spectral density function | equation (39)) into equation 
(GO) yields 



fr 
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u 1 
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IT 

a -ii) ! 



_ (l j <>') ? ( 1 - j 0 


I / 


(C7) 


Furthermore, substituting equation ((17) into equation ((11) yields 


2gM ] 

•* [ 

i -d d.-.j 

•f (1 -iOll 

u' 1 

o 

i - /* l 


1 Jl 
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The nondimcnsional exceedance rate \c 

o 

process is plotted in Figure r» as a function of 
ranging from 0.90 to 0.999. 


of the envelope of the g jitter 

A «r for various values of d 
g 


The values of t > indicated in Figure f> were selected such that the 

exceedance curve of tin* cnvcln|>c process remains below the exceedance curve 

of thc^ g jitter process (Fig# 2) Ibr the range of values of g <r and A *o 

g g 

indicated in Figures 2 and r>, respectively. This requirement results from the 
fact that the exceedance rate* of the envelope process should be less than the 
exceedance rate of tin* actual process at any prescribed level g a A o • 


g 


g 


The fact that the theory predicts exceedance rates of the envelope process which 
arc greater than the' exceedance rates of the actual process is a result of the 
fact that (filiation ((H) is strictly valid only in the limit 1 (i#e#, mono- 

chromatic g jitter)# Kquaticn ( (is ) is valid only as an approximation for d / 1 

t ♦ 

for A <r ■ A f <t , where A' • A . The quantity A is defined as 
g g 


Figure r>, Nonditnonsional expected rate 2»M -u* ^ of the envelope A(t) of tho 
g jitter process exceeding lovol A ^ o with positive slope for g > 0 ami 

negative slope for g < 0. {The expected rate of the envelope process 

exceeding level A o with positive or negative slope is 2M . 

g R 

The dashed lines correspond to tho condition 
A o > A’ , equation (hi)) . | 


Ji . 


1 - o' 


*(i - or 


and is that value of A e such that M N « 1 . We shall assume that A* = 

+ R R R + 

A 10. The Table contains a listing of the quantity A e for various values of 

R 

The dashed lines shown in Figure f> indicate those portions of the curves for 
which A/e -* A’ o . 







TABLE. THE QUANTITY A + /tr FOR VARIOUS VALUES OF 

g 


a 

+ , 

A /a 
g 

0 

0.798 

0.1 

0.934 

0.2 

1.111 

0.3 

1.344 

0.4 

1.661 

0.5 

2.111 

0.6 

2.793 

0.7 

3.936 

0.8 

6.232 

0.9 

13. 135 

0.95 

26. 95 

0.99 

137.5 

0.999 

15544.4 


D. Risk of Body Force Envelope Exceeding a Critical Value 

The risk R that the envelope of the g process will exceed level A 

from below for g > 0 or level -A from above for g < 0 at least once during 

c 

the orbital experiment time T is now given by 


-M T 

R = 1 - e K 


(70) 



where we have set X - 2M^ . Elimination of the g jitter envelope exceedance 
rate between equations (65) and (68) yields 


-fn(l - It) = 



/2\ Vz 1 - />' 3 
W (3(1- f.i) 


(n /?) 2 
1 


1/z -A*/2 <t 2 
e C g 
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(71) 












where we have substituted for ^/a with equation (67), Figures 6 through 

10 contain plots of A /a as a function of a.' T for various values of R and /3 
eg o 

ranging from 0.90 to 0.999. Again, the strong dependence of the statistics of 
the envelope process on ft is reflected in ihi e figures. It should be noted that 
because of the restrictions on the value of A for which the envelope exceedance 
analysis is valid, the parameter (i was restricted to values in the interval 
0.90 < < 1 for the construction of Figures 6 through 10. 



Figure 6. 


The quantity A /a 
c g 


as a function of a' T and R for / i = 0. 90. 
o 
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IV. SPACECRAFT ROTATION RATE EXCEEDANCE STATISTICS 

In this section we determine the exceedance statistics of each component 
of the 0 process referenced to the principal axes of the spacecraft, i.e. , the 
(x,y,z) frame discussed in_Paragraph II. C. 3. Our analysis will be concerned 
with a single component of fl ; consequently, we will dispense with subscripts on 
Q to denote components. We will use the symbol fl to denote a component of 
Q and cr^ to denote the standard deviation of fi . The developments in the 

subsequent sections are similar to those in the previous sections for the g 
process. However, the exceedance statistics of the envelope of the S2 process 
are not included because the available theory in the literature is valid only for 
random processes characterized by autospectral density functions which are 
symmetric about frequency u> m on the half-interval 0 < u> < ® ; i.e. , 

- u> m ) = 0 (oj + to^). The S2 process herein does not satisfy this condition. 




(73) 


where we have substituted equation (1G) into the integral to obtain the result 

indicated on the right side of equation (73) . The result indicated corresponds 

1 / 

to the geometric mean frequency (co to ) 2 with units of radians per second. 
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The ratios of the zero-crossing rates of , T, and g are given by 


n t.o N g.O 


a(i -y 




(74) 


N 

N 


T,0 

g»0 


= 1 


(75) 


Figure 11 provides a plot of N ft /N A as a function of (i according to equa- 

tion (74) which shows that N < N for all relevant values of (3. The 

St»u g*u 

reason for this result can be traced to the fact that g(t) is characterized by 
a flat spectrum over the domain Wj < u’ , while the spectrum of n(t) 

decreases as to -2 over the same frequency domain. This means that as |a>| 



Figure 11. The ratio of the zero-crossing rate of S2(t) to the 
zero-crossing rate of g(t) as a function of the 
bandwidth parameter /?. 
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increases from w, to u> the Fourier components of the g process will provide 
1 o 

increasingly larger contributions to the zero-crossing rate of g(t) , while the 
corresponding contribution to the zero-crossing rate of Sl(t) from each Fourier 
component of the SI process will be the same over the bandwidth dw for any 
frequency in the interval < I a? | < u.' o . 

Combining equations (72) and (73) yields the nondimensional exceedance 

rate 


2jiN -Sl 2 /2or 2 

—2 = o 4 c 0 

CO 

o 


(7G) 


Figure 12 contains a plot of 2?rN_ /a' us u function of SI / a for p’ = 0.01 and 1. 

SI o Si 



Figure 12. Nondimensional expected rate 2 ttN^ 4*?^ of exceeding 

the nondimensional level |Si | with positive slope for si > 0 

or negative slope for S2 < 0. (The expected rate of exceeding 
level |S1 I /o with positive or negative slope is 2N .) 

C Si Si 



B. Risk of Rotation Rate Exceeding a Critical Value 

Following the developments in Paragraph m.n, wo seek to determine 
the risk that the In I process will exceed a critical value In I for a given 

orbital experiment duration time T. We hypothesize that the number of exceed- 
ances of In | above level |n | from below is a Poisson process. Thus, the 

c 

risk R that the In I process will exceed the critical value In I < at least once 
during an experiment of duration time T is 




This formula permits the calculation of a critical rotation rate |n I as a func- 
tion of risk U of the quantity In I exceeding level In I at least once during an 

experiment of duration time T and the n process spectral density parameters 

a' and u 1 . Figures 13 and 11 contain plots of In I V. as a function of Tu’ 
o n o 

for various values of risk for /> 0.01 and 0.990, respectively. 


V. CONCLUDING COMMENTS 

The previous sections provide a simple stochastic moiiel of spacecraft 
rotation and induced gravity. To develop the model, it was assumed that the 
components of the net applied torque vector are stochastically independent 
Gaussian processes. Validation of this model must await the results of 
statistical analyses of accelerometer and rate gyro data acquired from past 










spaceflight missions. The authors of this report are currently analyzing thruster 
rate gyro and accelerometer da( i acquired on the Apollo-Soyuz mission. However, 
in the interim time period the proposed model can be used in orbital experiment 
definition studies. If it is found that the components of induced gravity and 
vehicle rotation are non-Gaussian processes, then the calculation of risk values 
associated with assigned critical values of vehicle rotation and induced gravity 
could prove to be an extremely complex task for future spaceflight missions. 

To apply the model to obtain estimates of rotation rate and induced gravity 
spectra and risk values associated with exceeding critical values of rotation rate 
and induced gravity, estimates of (i» u) q , o , and are required. A range of 

values for each parameter should be used to obtain a ’’feel" for the effects of 
g and ft on an experiment. It should be remembered that the standard devla 
tions of the components of the rotation and induced gravity vectors are related. 

In fact, the standard deviations and cross-variances of the components of g 
are derivable from the standard deviations of the components of ft via equations 
(18) and (28) through (34) upon specification of the vehicle principal moments 
of inertia and experiment location relative to the vehicle center of mass. 

It should also be remembered that the model described is valid for a 
particular imposed torque process; i.e. , the values of /}, u' q , and the standard 

deviations take on fixed values. However, these quantities can vary 
in time during a mission. To obtain exceedance statistics of g and ft and associ- 
ated risks of exceeding critical values of S2 and g for this case, uk* j >int probability 
density functions p.(/j, w , o 0 ) and p 0 (/i, u , a J for the mission are required. Con- 

sidering the relationships between the standard deviations of the components fi and g, 
these iunctions should be derivable from one another. The expected exceedance 
:'ates of ft and g for a total mission can l>e obtained from the following 
integrals: 


00 00 1 



N (SI) 
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CD 00 J 
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where N in the integrands of equations (79) and (80) are given by equations 
( 56) and ( 72) , The calculation of risks of exceeding critical values of g and Q 
can be performed by invoking the Poisson model used in the previous sections. 

A similar analysis can ac applied to the expected exceedance rate of the envelope 
of the g process. In this model for the exceedance rates of g and SI , the torque 
vector process is assumed to be piecewise statistically stationary in time. 
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